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Abstract — Opportunistic routing is a multi-hop routing scheme 
which allows for selection of the best immediately available relay. 
In blind opportunistic routing protocols, where transmitters 
blindly broadcast without knowledge of the surrounding nodes, 
two fundamental design parameters are the node transmission 
probability and the transmission spectral efficiency. In this 
paper these parameters are selected to maximize end-to-end 
performance, characterized by the product of transmitter density, 
hop distance and rate. Due to the intractability of the problem 
as stated, an approximation function is examined which proves 
reasonably accurate. Our results show how the above design pa- 
rameters should be selected based on inherent system parameters 
such as the path loss exponent and the noise level. 

I. Introduction 

Ad hoc networks operate on the basis of multi-hop rout- 
ing, which allows information to be communicated across 
the network over a series of hops. End-to-end performance 
depends critically on the quality of the multi-hop routes used, 
but choosing good routes in a dynamic network (e.g. one 
where nodes are moving quickly) is a particularly difficult 
task. Opportunistic routing (e.g fl] - ||6l) is well suited to such 
dynamic settings, as it can be performed with low overhead 
while also exploiting spatial diversity (in terms of fading and 
topology). 

In this paper we focus on a purely opportunistic routing 
protocol, in which nodes broadcast packets in a completely 
blind manner (i.e. without any knowledge of surrounding 
nodes). More specifically, the protocol is described as follows: 

1) In each time slot, each node randomly decides, with 
some fixed probability, to either transmit or receive. 

2) A transmitting node will blindly broadcast a packet at 
some fixed spectral efficiency, whereas receivers listen 
for any transmission. 

3) All receiver nodes that successfully decode the packet 
will send an ACK back to the transmitting node during 
an acknowledgement period. The ACK packets contain 
absolute geographic information about the RX location. 

4) The transmitter selects the successful receiver which 
offers the most progress towards the packet's final des- 
tination and sends a message to that node electing it to 
become the next forwarder. 

This protocol allows for the instantaneous choice of the 
best next hop, at the cost of an acknowledgement period and 
with the requirement for geographical information. Within this 



algorithm, the key adjustable parameters are the transmission 
probability and the transmission spectral efficiency. 

The transmission probability, denoted by p, determines the 
proportion of transmitters to receivers in each slot. When p 
is large, there is a large amount of interference between the 
nodes. This interference will cause only receivers close to 
transmitters to be potential forwarders. When the transmission 
probability is low, there are fewer simultaneous transmissions 
and thus less interference and more available relays. Thus 
longer hops are possible, but in fewer numbers. As a result, the 
trade-off with transmission probability is essentially between 
many simultaneous short hops, or fewer long hops. 

When the spectral efficiency is high, a large signal-to- 
interference ratio (SIR) is required to decode. Thus, only relays 
which are close to the transmitter are likely to decode. On the 
other hand, a lower spectral efficiency allows nodes that are 
farther away to decode. Therefore, the trade-off with spectral 
efficiency is between shorter hops at higher data rate or longer 
hops at a lower data rate. 

The objective of this paper is determining the transmission 
probability and spectral efficiency that optimally balance these 
trade-offs. We study blind opportunistic routing in a spatial 
model, as in In this model, end-to-end performance is 
characterized by the forward-rate-density, which is the prod- 
uct of average transmitter density, average hop distance and 
transmission rate (spectral efficiency). This is a function of 
measured network properties (e.g. path-loss exponent) and 
adjustable parameters (spectral efficiency and transmission 
probability). Our objective is to find the transmission prob- 
ability and transmission spectral efficiency that maximize this 
metric. Studying this quantity directly is intractable, however 
we derive a reasonably accurate approximation and reinforce 
the conclusions with Monte Carlo simulations. 

II. Preliminaries 

A. Network Model 

Consider an infinite set of transmit/receive nodes $ dis- 
tributed according to a homogeneous 2-D Poisson point 
process (PPP) with density A m^^. We consider a slotted 
transmission scheme. In each slot a node elects to become a 
transmitter with probability p, independent across users and 
slots. The set of locations of the transmitters (TXl, TX2, 
..) denoted <&* and the set of locations of the receivers 
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Fig. 1 . Example network: TX node 1 has a message to be sent along dashed 
line. Among successful message receivers, RX node 2 offers the most forward 
progress and will forward the packet in the next slot. 



then form independent PPP's of intensity Xp and A(l — p) 
respectively. 

B. Channel Model 

We consider a path-loss model with exponent a > 2 and 
Rayleigh fading coefficients ft,, j from TX i to RX j. Denoting 
the signal transmitted by TXi (located at Xi) as Ui, the 
received signal of RXj (at location Xj) is given by: 



Y, 



Xi — X 



l-a/2 



Ui + Z. 



(1) 



For the time being, consider the network is interference 
limited, so thermal noise (z) is negligible. The case of a 
channel with thermal noise is considered in Sec HV] 

Assuming the transmit powers are all equal, the SIR (signal 
to interference ratio) from TX i to RX j is defined: 
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\Xk — Xi 



(2) 



We assume that all users transmit at rate equal to R{/3) = 
log2(l + /?); thus a communication is successful if and only 
if the received SIR is larger than a threshold value /3. 

III. Progress-Rate Product 

Although we are interested in the network-wide perfor- 
mance of the protocol, by the homogeneity of the PPP this 
is statistically equivalent to the performance of any transmit- 
ter. Thus, without loss of generality we can focus on TXO 
located at the origin. From the perspective of TXO, all the 
other TX nodes form a homogeneous PPP with intensity Xp 
(by Slivnyaks Theorem |7|, the distribution is unaffected by 
conditioning on the presence of TXO). The transmitters and 
receivers are ordered by their distance to the origin i.e. TXI 
is closest to the origin, TX2 the next, etc. For the sake of 
simplicity, the SIR from TXO to RXj is denoted Sj and is 
given by Sj = Sqj. 



A. Forward progress density 

We assume that the transmitter located at the origin has 
a packet to be sent to a receiver a very large distance away 
(e.g., fig. [T]). Specifically, we assume the final destination is 
located an infinite distance away along the x-axis. Thus, the 
relay offering the most progress is the one with the maximum 
x-coordinate amongst the set of receivers that successfully 
decode TXO's packet. The corresponding progress, for a given 
network realization, is: 



max 



1(5, >/3|$*)(|X,|.cos (eix,)) 



(3) 



where 1(5', > /3|<i>*) is the indicator function that the SIR at 
RXj (at point Xj) relative to TXO (at the origin) is greater than 
the threshold. The expectation of this function with respect to 
the node process and fading is d{p,j3). 



dip, (3) 



(4) 



This quantity is the expected progress (towards the destination) 
made in each hop. The forward progress density f{p,l3) then 
is: 

f{p,(3)^d{p,(3)Xp. (5) 

Finally, the forward progress-rate-density function K{p,/3) is 

K{p,f3) ^ (Ap)•E[I?(<i>^<^)^/3)]•log2(l + /3) 

= /(p,/3)-i?(/3). (6) 

Assuming each node in the network wishes to send data to 
another node a distance L away, at most p ■ d{p, (3) /L packets 
can originate at each node in each slot UJ. Translating to 
bits/sec/Hz, this means that the maximum end-to-end data rate 
(per node) is 



p-dip,l3)-\og^{l + P)/L [bps/Hz] 



(7) 



Thus, K{p, P) is directly proportional to end-to-end rate. 

The quantity K{p, (3), which is deterministic, is a function 
of transmission probability p and SIR threshold (3. We seek 
the optimal values of /3 and p which maximize this product: 



{p*,l3*) = argmax7^(p,/3). 

P,I3 



(8) 



A closed-form expression for K{p, (3) cannot be found in 
general, so we rely on an approximation which is developed 
in the following section. 

B. SIR Cell-Based Approximation 

The basis of our approximation to K{p, f3) is the concept of 
an SIR cell fS], ID. In the absence of fading, any RX within a 
particular region (i.e., the SIR cell) around a TX can decode, 
and the best forwarder is the RX in this region with the largest 
a;-coordinate. Therefore, the forward progress is determined by 
the SIR cell - which is completely determined by the interferer 
locations and thus is a function only of the interferer process 
- and by the receiver locations (i.e., $''). To reach a tractable 
approximation for K{p, (3), we remove the randomness in the 



SIR cell and assume it is a deterministic region that is a 
function of p and (3, while retaining the randomness in <!>''. 

We consider the average size of the SIR cell. This quantity 
is clearly defined in the absence of fading. On the other hand, 
the SIR cell is not well defined with fading because the SIR 
depends on independent fading RV's (specific to each receiver 
location) in addition to the interferer locations. However, we 
can derive a quantity analogous to the SIR cell area, denoted 
by vq, by integrating the point- wise success probability over 
all space. From [9 |, the success probability for a RX a distance 
y from the TX is: 



Po(y) 



Thus, vq is given by: 
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and r(z) is the gamma function r(2:) = t^^^e^^dz. 

We will then assume the SIR cell has area equal to its mean. 
We now approximate the SIR cell by a square centered on TXO 
with side ^JvQ. In other words, we assume that all RX within 
this square are able to decode (ref. fig. Vf^ is the positive 
half of that square with area vo/2. Under these assumptions, 
the expected forward progress is: 
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where #(•) gives the cardinality of a set and 

\-{l-p)-vo _ (1-p) G[a) 
^ 2 p 2^2/a- 

Note that the value c is the expected number of receivers in 
the half-square V(^; it is also equal to half the average number 
of successful outgoing transmissions from a TX as in [|4j. 
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Fig. 2. SIR cell is assumed to be a square Vq of side length ^/vq. The 
probability that the forward progress is x is the probability there are no RX 
in the dark shaded area. 



Thus, our approximation to the progress-rate-density is: 



^(p,/3) = 4,(p,/3)(Ap) 



(14) 



K,,{p,P) = f,,ip,(3)\og^il+/3). (15) 



The optimal values then are the solution to the equation 

0. (16) 



Although there is no closed form solution to this equation, 
the approximation nonetheless yields valuable insight. By rear- 
ranging the terms in ( fT2] i. we can interpret the forward progress 
approximation as the product of the maximum possible relay 
distance (i.e., half the side-length of the SIR cell) and a 
fractional term that captures how much of the maximum is 
attained by the best receiver 
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max. relay dist. ^^^^ „f 

The value -yo is decreasing in p: as p increases, the amount 
of interference between transmitters is increasing, giving each 
transmitter a smaller expected SIR cell size. The value of c is 
also decreasing in p as a result of decreased wq and smaller 
receiver density. 

The value Kaq{p, /3) can be written: 



Ksqip,/3) 



log2(l + /3)VApG(a)/3 
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(18) 



which lends the intuition that the end-to-end performance is 
increasing with p until the 1 — ^—^ — fractional term becomes 
too small. When the value of c (expected potential forwarders) 
is greater than 3, this fractional term is greater than 2/3. When 
p is large enough to cause c to go below 3, the performance 
drops off sharply. 

IV. Numerical Results 

For numerical results, Monte-Carlo simulations were per- 
formed. Locations of receivers and transmitters were realized 
using the properties of the PPP (c.f. |10|). SIR was calculated 
at each receiver and the largest progress offered was averaged. 



A. Accuracy, Path-loss and Fading 

Fig. [3] shows the accuracy of the approximation of the 
spatial density of progress (|5]), QHl for p and (3. While the 
approximation is not absolutely accurate, it does correctly 
capture the dependence on p and (3. 

Fig. |4] shows the best p or /3 maximizing the progress- 
rate-density (fTST i for a fixed /3 or p, respectively. Again, the 
approximation helps to answer the essential question of what 
parameter choices are optimal. 

In Fig. |5] the jointly optimal {p,(3) are plotted for different 
values of the path loss exponent. Note that the optimal value 
of p is essentially the same for a = 3 and 4, but that 
the optimizing SIR threshold increases with the path loss 
exponent. This increase is similar to im . where the optimizing 
SIR threshold is found for a slightly different, but related, 
problem. 

B. Sensitivity to parameter variation 

Fig. |6] shows the normalized approximation of Eqn. ( fTSl) for 
path loss a — 3 defined by 



-^sq,iiorm(P5 /5) — 



(p,/3) 



(19) 



This figure shows how the performance decreases with differ- 
ent values of p and /3. Observe that there exist a wide range of 
(p, /3) pairs that are near optimal (within 90%). If there exists 
a strong reason to choose a particular /3 (e.g., limited code 
rates & modulations), or a particular p (e.g., power cycling, 
energy saving from different TX and RX powers), end-to-end 
performance does not suffer as long as the other parameter is 
appropriately chosen (to figs. H] |6]l 

C. Robustness to Noise 

In order to simulate the effect of noise on the forward 
progress, we modify the SIR ^ to include noise power a^. 



Si J — 



\^\Xk-X, 



(20) 



The metric for comparison is the average SNR to a node 
located at the average nearest-neighbor distance {(Inn)- 



SNRnn - ^iml ^ (21) 

Fig. |7] shows the maximum forward-rate-density is decreasing 
with increasing channel noise. Also, the number of potential 
forwarders is decreasing with increased channel noise. 

Fig. [U shows the effects of channel noise power on optimal 
p and [3. The value of p* is increasing and the value of (3* 
is decreasing with higher noise floor The maximum possible 
spectral efficiency is expected to drop with increasing noise 
power As explored in IIV-BI when /3 is fixed to some value 
(in this case as a result of the noise floor), the transmission 
probability can increase to compensate. 

' Simulation results are similar for a non-fading environment. 
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Fig. 3. Spatial density of progress: approximation (thick curves) and 
simulation for A = 1, a = 3, varying SIR threshold /3. 
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Fig. 4. Optimal p values giving best progress-rate-density for fixed values of 
SIR 13 (top). Optimal /3 values giving maximum progress-rate for fixed values 
p (bottom). 
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Fig. 5. Optimal choice of p and /3 giving maximum progress-rate for fixed 
a path-loss values. 
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Fig. 7. Simulated best progress-rate-density achievable for varying SNR to 
node at mean nearest-neighbor distance (A7r)°/^/(T'^ 



We studied the case of opportunistic routing in an ad hoc 
wireless network with the goal of maximizing the product of 
forward progress density and the rate of data transmission. We 
developed an approximation using the concept of SIR cells 
and found for a — 3 the optimal spectral efficiency is near 
1.3 bps/Hz and the optimal probability of transmission is near 
0.06. The results give the intuition that having a low transmis- 
sion probability ensures a low amount of interference and a 
relatively high SIR threshold ensures only local transmissions 
are received. 
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